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Abstract 

We use the Matrix formalism to investigate what happens to strings 
above the Hagedorn temperaure. We show that it is not a hmiting 
temperature but a temperature at which the continuum string picture 
breaks down. We study a collection of D-O-branes arranged to form 
a string having N units of light cone momentum. We find that at high 
temperatures the favoured phase is one where the string world sheet 
has disappeared and the low energy degrees of freedom consists of 
A^^ massless particles ( "gluons" ) . The nature of the transition is very 
similar to the deconfinement transition in large-N Yang Mills theories. 
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1 Introduction 



A feature of string theory that has not been fully understood is the exis- 
tence of a "limiting" temperature, often called the "Hagedorn" temperature 
H, 0, 0|. The reason for a limiting temperature is the exponential growth 
with energy of the number of states in a relativistic string which makes the 
canonical partition function diverge at this temperature. A lot of work was 
done in exploring various aspects of this phenomenon [||§]-[13]]. The physical 
picture that has emerged is roughly as follows . Consider a gas of weakly 
interacting strings. As the Hagedorn temperature is reached it is entropically 
favoured for most of the energy to go into one big string. As a result the 
temperature of the string gas (measured by the average kinetic energy of the 
strings) does not increase - hence the idea of a limiting temperature. The 
problem is subtle because, as noted above, the canonical ensemble does not 
make sense near the Hagedorn temperature. The partition function diverges, 
the specific heat is negative and a microcanonical description has to be used 
to arrive at this picture. 

What is not clear from this analysis is whether the Hagedorn tempera- 
ture is really a limiting temperature or whether it is to be interpreted as a 
phase transition temperature. One motivation for interpreting it as a phase 
transition comes from the belief that large-N Yang-Mills theory has a string 
description, where the string coupling constant, Qs, can be identified with 
This theory is believed to have a transition (for any A^) from a low 
temperature confining phase where the lightest degrees of freedom are colour 
singlets, to a high temperature deconfined phase where the light degrees of 
freedom are gluons (and quarks). It is only natural to identify the Hagedorn 
temperature with the temperature of the deconfinement transition. 

Evidence for the existence of two distinct phases in Yang-Mills theory, 
in the large A^ limit, was given in WM. It was shown that in the high tem- 



perature phase asymptotic freedom requires that there be a term of 0(A^^) 
in the free energy which is absent in the low temperature phase where the 
degrees of freedom are colour singlets. It was also noted there that a possible 
explanation for the different form of the expression for the free energy is the 
non-analyticity caused by an exponential (in energy) increase in the number 
of states in the low temperature phase leading to a Hagedorn temperature. 

A further step in this direction was made in [^. There is an interpre- 
tation of the Hagedorn transition as a Kosterlitz-Thouless transition on the 
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string world sheet induced by a condensation of vortices |T0[. The vortices 
correspond to strings wound around a compact Euchdean time coordinate. 
The Kosterhtz-Thouless point is where this winding mode becomes mass- 
less. Using this picture it was argued in [Q, that the expression for the 
free energy {—InZ) acquires an additional term of O(p-) = 0{N'^) after the 
phase transition. Using the results of |n| this made the connection with the 
deconfinement transition stronger. Furthermore this extra term is a temper- 
ature dependent genus zero contribution, which cannot ordinarily be there if 
a continuum string world sheet exists. Thus the world sheet must somehow 
disappear (this can happen due to a condensation of vortices) and this would 
correspond in the gauge theory to a liberation of gluons. Thus a consistent 
physical interpretation of the Hagedorn transition, based albeit on indirect 
arguments and analogies, can be made. 

Issues such as the phase above the Hagedorn transition where the string 
supposedly disintegrates cannot really be discussed using string perturbation 
theory. However, in the last year or two a formalism that purports to treat 
string theory in a non-perturbative fashion has been proposed P^l 53] and 



various issues have been studied [[Q" IHI- Finite temperature aspects of 



this formalism have also been discussed It is appropriate therefore to 
investigate the Hagedorn transition once again in the light of new insights 
that have been obtained using this formalism. This is the subject of this 
paper. We will show that in the matrix model of N D-O-branes, at low 
temperatures, there is a string phase where the D-O-branes are spread out to 
form a membrane wound around a compact direction. (This picture has also 
been used in ||3^ to discuss other aspects of the Hagedorn transition.) In this 
phase the ground state must be that of a HA string with a specified value of 
string tension. The light modes are those of a string. At high temperatures 
this string ceases to exist. The light degrees of freedom are massless modes. 
These are just fluctuations about the origin of the A^^ matrix elements. This 
phase transition can be identified with the Hagedorn transition for the HA 
string constructed out of the A^ D-O-branes. [] 

The matrix model calculation of free energy in the low temperature phase 
is suspect because we do not know the exact zero-temperature expression. 

^This issue has also been addressed in the "string-bits" formaUsm, with similar con- 
clusions 1^ . The essential idea being that the string is a polymer of "string-bits" , and at 
high temperature dissociates into its constituents. 
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In particular one knows that there must exist a threshold bound state for 
any N (this has been proved for N = 2[^, ^) and this is not seen in 
perturbation theory. Our calculation of the (temperature dependent part of) 
the effective potential which is only a one-loop approximation is therefore 
not good enough in the low temperature phase. However, even though non- 
perturbative effects are expected to change the number of light degrees of 
freedom, one can argue that the way it does this is by reducing the number 
of flat directions and this can only reduce the number of light degrees of 
freedom. Thus hopefully we have overestimated the magnitude of the free 
energy and therefore the existence of the phase transition can still be argued. 

In matrix theory a membrane is described in the limit of being large 
by the following configuration 



X' = Up; = Uq (1) 
Here p, q are matrices that satisfy 

[p. = ^ (2) 

One can consider a representation (as in [^]) where the eigenvalues range 
from to 27r in steps. (^ represents a membrane in the ij plane, of sides 
2'kU and 2txU . If i,j happen to be compact dimensions of radius U'\ then 
(|I]) describes a wrapped membrane, which classically describes a stable BPS 
configuration. We can also consider a background configuration of the form 



m 



X' = Vp (3) 

which describes a membrane wrapped around X* with the other edge free i.e. 
a free string, not necessarily wrapped. If we let the X^ be periodic functions 
of q of the form e*"^'', then we get a closed string. Let us construct this string 
action in the matrix model. 

The matrix model for M-theory |^ is essentially a D-O-brane action. 



reinterpreted as a DLCQ M-theory where the "eleventh" dimension, of 



radius , is a light cone direction. This is related to i^n of the M- 



theory that we would like to study by a very large boost. The component 

jv 



of momentum is (where N is the number of D-O-branes), and is like a 



^Explicit representations can be easily found in the literature, e.g. H, Esl 
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ten dimensional mass. When P"*" is large (as it is in an infinite momentum 
frame) the action is that of a non relativistic particle with mass Thus 

the evolution operator = *^p+ is the Hamiltonian obtained from the 
D-O-brane action. The bosonic part of the action for N D-O-branes is 

We have explicitly written factors of Ig, to make the action dimensionless. 
Our conventions are as follows: Qs is the (HA) string coupling constant and 
we have defined the inverse string tension to be 27ia' with a' = Is ■ The 
parameters of DLCQ M-theory are the radius R~ and the eleven dimensional 
Planck length 1^^. 1^^ is defined by the membrane tension which we have taken 

to be — • This fixes = (ttt)'^ and also a' = £_ . These relations 

also imply that gsh = R ■ Thus the kinetic term is essentially "|mf^" with 
m = The coefficient of the potential term can be seen to be ^,11 . The 

normalization of the potential term is chosen so that the classical mass of a 
D-2 brane comes out right. The parameters of M-theory are related to those 
of this DLCQ matrix theory by the scaling relations ||] -^jtt^ = jlf > 

L9 _ L9 

W ~ /II ■ 

> '■p 

Let us assume that is a compact dimension of radius L^, assumed to 
be small. We will take this as the "eleventh dimension" . When a membrane 

is wrapped around we get a string with inverse tension f3' = ^9 , and 

string coupling g^^^, = 

The 1 + 1 dimensional action that we obtain is that of 9+1 dimensional 
Supersymmetric U{N) Yang-Mills theory reduced to 1+1 dimension. 



+ {F„,y- + ^,{lX',.P]r] (5) 



This supersymmetric action has also been studied in a different context by 
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Dx = dx + iA^ is the covariant derivative in a direction X^*, which is 
T-dual to X^, and is of radius Lg. x is thus a coordinate along a D-l-brane 
wound around X^*. If X^„ is a matrix describing the lattice of D-O-branes, 



(m,n = 0...M) then following we have 



1 r 

= ^ E e-^^ol (6) 

Note that Xqo is the original D-O-brane matrix of the uncompactified theory. 
For us Xqq will be given by (0). Thus Dx is given by 

dx®I + I® 4^9, (7) 
a iV 

The derivatives act on eigenf unctions: 

e ^e''^Pe''''' (8) 
The eigenvalues are "^^tl" . Thus the effective radius is XLg. These are the 



long strings described in qO|, |2^, We can replace D^; by da and 

the eigenfunctions by e'^T^ . Here a has a range of — 2t[NL . Using 
Lg = ^ = we can rewrite the action (||): 

2wN0' 

dt [ da{{dtX' f - {daX' f + ...] (9) 
Jo 



Noting that is P"^, we recognize the light cone gauge string action with 
a string tension of (3' . Turning on Fqq corresponds to addition of D-strings 
[^ . The commutator terms are zero if we restrict the matrices X* to be 
X*(a;, i.e. without any p dependence, p dependence corresponds to 
fluctuations in the matrix model that are not string-like. If we knew the 
fully non-perturbative effective action presumably it would be manifest that 
these fluctuations are massive. 

If one calculates the partition function using the above string action one 
will encounter the usual divergence as the Hagedorn temperature is reached. 
However we see that the divergence is an artifact. At the Hagedorn temper- 
ature energy in the string is overwhelmingly found in the highest mode that 
can be excited with that energy. In the matrix model there is a cutoff on 
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the mode number which is in fact NM. Thus in the partition function there 
is a natural cutoff in the energy integral above which the discrete nature of 
the string becomes important and at this point one has to go back to the 
original matrix model. Thus we see that the Hagedorn temperature is not a 
limiting temperature. Of course eventually one has to let NM — > oo (in fact 
one should sum over all values of A^). Presumably some sort of scaling has 
to be done to extract finite physical quantities, as has been discussed in, for 

e 



example, [^, ^ We should also keep in mind that j3' = P'i-pj 



Thus — > oo as /p"*^ — > 0, in order to get a finite 

Let us try to understand the nature of the phase transition. We have 
shown above the emergence of a string in the matrix model when the eigen- 
values of a particular X* are spread over a finite range. L can thus be thought 
of as an order parameter characterizing the two phases of the model. When 
L is non-zero we have a configuration where the eigenvalues of X* are uni- 
formly distributed from to 2ttL and it describes N D-O-branes arranged to 
give a string, whereas if L is zero we get a cluster of N D-O-branes located 
at the origin and no string. Given this, there is a simple way to understand 
the nature of the phase transition. The issue of whether a uniform distribu- 
tion of eigenvalues is the preferred state or not goes back to the issue of the 



breaking of "f/(l) symmetry in the reduced Eguchi-Kawai models [[|T5| 



A phase transition of this type as a function of the coupling constant 



in large-N lattice gauge theory was also noted in pOl pTl. In non supersym- 



metric models even in the one loop approximation one can see that there is 
a logarithmic attractive potential between the eigenvalues and they tend to 
cluster, thus breaking the U{lY symmetry. In supersymmetric models on 
the other hand, Bose- Fermi cancellation allows a uniform distribution [P^] . 
Thus if one assumes that supersymmetry somehow prevents the clustering 
of eigenvalues, then one is led immediately to conclude that at finite tem- 
perature, supersymmetry being broken, the eigenvalues will start clustering. 
This then is the basic mechanism of the Hagedorn transition. 

Of course the above argument does not explain why there should be a 
finite transition temperature - supersymmetry is broken for arbitrarily low 
temperature. This clearly has to do with the uniqueness of the assumed 
ground state. Somehow the flat directions should be lifted and a non-trivial 
potential generated. This will ensure that only at a finite temperature will 
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the temperature effects dominate and change the nature of the equihbrium 
(themodynamic) state. Thus one needs to know the exact potential to make 
this argument concrete. 

We do not know how to compute the exact potential but it is sufficient 
for our purposes that the configuration should have non-zero overlap with 
the (unique) ground state and therefore lies in the "bound " region of what- 
ever non-trivial potential for L the exact theory has. The argument for this 
is as follows: We know that for D-O-branes there is exactly one threshold 
bound state. We also know that this is the (super)graviton with units 
of momentum along R^^. From the viewpoint of the 10 dimensional theory 
where is the eleventh direction and X^^ is the light-cone direction, this 
is the graviton with A^-units of light-cone (p^^) momentum. But this gravi- 
ton is precisely the ground state of the light-cone string described by (|^). 
Therefore this configuration must have a non-zero overlap with the unique 
non-degenerate ground state. Thus if L denotes the spread of eigenvalues, 
which is the size of the bound state, the ground state wavefunction denoted 
by ip{L) is non zero for a range of values of L and therefore < L'^\il) >^ 0. 
Here the ket \U > denotes the state corresponding to the classical configura- 
tion given by (^. The above considerations thus explain why the Hagedorn 
transition happens at a finite temperature. 

We now turn to a computation of the temperature dependent free energy. 
It will be a one loop computation about two different backgrounds, one where 
L is non-zero and one where L = 0. We will assume that one dimension is 
compact, say , with radius L^. We will use the standard trick of compact- 
ifying the time coordinate t in the matrix model to a range of (3 (not to be 
confused with the inverse tension (5'). Of course t is conjugate to P~ rather 
than H, but in the infinite momentum limit for the 10-dimensional theory 
it is the non-relativistic energy and therefore t is a non-relativistic approx- 
imation to the time coordinate. The calculation of InZ is straightforward 
(||23||. We have to keep in mind that fermions have anti-periodic boundary 
conditions in time. The bosonic contribution is, 

5Tr In E (10) 

The fermionic contribution is 

- 4Tr In E (11) 
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and the ghost contribution is 



- ITr In E (12) 

where E = Pq + P| + 6^. Here Po = j for bosons and ghosts and Pq = 12-^/^ 
for fermions. Pg = ^ for all the fields. R is some length (which could be 
L*^ or NL^* for instance), b is some small parameter needed to regularize 
divergences and which will normally be set to zero. Thus we obtain 

4iVE / — (e"^^+^^' - e"^^^+^^') (13) 







The factor comes from the degeneracy due to the factor e*''^ (where r has 
a range from to — 1) in the eigenfunctions 

Using a Poisson resummation on m one reduces this to 



p odd ' 

= -SAT E -e~'"«^" (14) 

p odd ,n P 

When (3 oo, terms with n ^ clearly go to zero. When n = one has to 
introduce the regulator b (which replaces in the above formula) and take the 
limit 6 — > after taking the zero temperature limit and again the expression 
vanishes. This is as it should be indicating that at zero temperature Bose- 
Fermi cancellation makes the one loop correction vanish. 

The zero mode term n = gives a logarithmic divergence due to the sum 
over p. On regularizing with b we get 8Nln bjS. For large ^ the rest of the 

contribution is —8Ne~^. R in the low temperature case is NL^*. Thus we 
get 

/3F = N{8ln 13b - 8e'^) + higher orders (15) 

The logarithmic term can be seen to be the classical (high temperature) term 
in the free energy of a harmonic oscillator (any temperature is high relative 

'''As mentioned earlier the factor N is an overestimate. The exact potential would break 
the remaining degeneracy and should make it 0(1). 
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to zero!).Q The second term is just the leading term in a low temperature 
expansion of the free energy of a bunch of harmonic oscillators which goes to 
zero as /5 — s> 0. Thus at low temperatures it behaves like a zero-dimensional 
point particle rather than a one-dimensional field theory. 
At high temperatures (small P) we get 

PF = N{8ln m + 8ln^- ^) + ... (16) 

K 2 

The leading behaviour is thus that of a classical one-dimensional field theory. 

We can consider two extreme cases: R = NL^*, which corresponds to a 
string configuration with the order parameter L being equal to L^*. In that 
case (0) becomes 

(3F = N{2ln m + lnj^^ ^) + ... (17) 

We can see the connection with string thermodynamics as follows: From 
the above result for the free energy, 

7-12 

<p- >= — ^ NL^*T^ 



and this gives E ^ \/^^T. Using S ^ NL^*T (and L^*R- = /?') we get 

S ^ ^/J3 E, characteristic of strings. 

In the case where all the D-O-branes are clustered at the origin and L = 0, 
we have to set R = L^*. There is one other difference with the previous case. 
The degeneracy of each level is A^^ rather than A^. This is because the 
eigenvalues of acting on the eigenfunctions given in (|^) do not depend on 
n,m (i.e. the eigenfunctions can have any p,q dependence). Thus the entire 
expression for the free energy is now multiplied by A^^. This gives 

(3F = N\2ln m + ln^- ^— ) + ... (18) 

^The zero mode corresponds to a free particle whose partition function (in one dimen- 
sion) goes as \/TmV where V is the (infinite) volume . However when we add the 6-term 
(a harmonic oscillator potential term) , it regulates this infrared divergence and contributes 

a factor of i / ^ (in place of V) where k is the "spring constant " , and the net contribution 



• T 
IS ^. 
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Notice that the leading high temperature piece is the same in both cases. 
Both are proportional to N"^. However the non- leading term is different. 
Clearly for small f3 when the logarithm becomes negative (]TB|) is more neg- 



ative than ([T7|). Hence the phase without the string is favoured at high 
temperatures. 

Let us summarize what we have learnt. The phase below the Hagedorn 
temperature corresponds to one where the D-O-branes are arranged to form 
a membrane which can be wrapped around a compact direction to give a 
string. At high temperatures the D-O-branes prefer to cluster at one point, 
thus the string disappears. We cannot estimate the temperature at which 
this happens without some knowledge of the exact potential as a function 
of L ^ If the height of the potential at L = is of order then the 



transition can be expected to take place at around a temperature of — ^ 

which is essentially the Hagedorn temperature. But there is no estimate of 
this in the literature as far as we know. 

The other thing we have learnt is that the nature of the phase transition 



is very similar to the deconfinement transition in large-N Yang-Mills. In [51 



it was shown that the deconfinement transition has exactly the same form. 
The two phases described there were the same: One where the eigenvalues 
are uniformly distributed and one where they are clustered at one value. 
Furthermore the fact that at high energies the free energy has a term of 0(A^^) 
and that it corresponds to massless particles in the adjoint representation 
("gluons") also makes the connection more concrete. 

The spacetime picture above this phase transition is still a little myste- 
rious. Does the cluster of D-O-branes at the origin still represent gravitons 



or perhaps a black hole? That is reminiscent of situations discussed in [52 



Also, the membrane of M-theory can wind around any compact dimension 
and give a stable string configuration. In particular it can wind around X^^. 
These strings are the ones that connect the D-O-branes. As the matrix model 
does not contain these string excitations, (indeed the scaling prescription of 



13, @ ensures this), we cannot discuss them. Nevertheless if we interpret 



the results of this paper to mean that the membrane configuration itself 
is unstable then there cannot be any strings left at all. It is also possible 
that there are other phase transitions at temperatures T ^ 0(1/ Ru), which 

^An attempt to calculate the potential is described in |39| 
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cannot be computed in this model. 
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